3rd Order Temporal Correlation Function of Pseudo-Thermal Light 
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This experiment reports a nontrivial third-order temporal correlation of chaotic-thermal light in 
which the randomly radiated thermal light is observed to have a 6-times greater chance of being 
captured by three individual photodetectors simultaneously than that of being captured by three 
photodetectors at different times (separated by the coherent time of pseudo-thermal light), indicating 
a "three-photon bunching" effect. The nontrivial correlation of thermal light is the result of multi- 
photon interference. 
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PACS numbers: PACS Number: 03.65.Bz, 42.50.Dv 

In 1956, Hanbruy-Brown and Twiss (HBT) discovered 
a surprising photon bunching effect of thermal light [![. 
In a spatial HBT interferometer, the spatially, randomly 
distributed thermal light was observed to have a two- 
times greater chance of being captured by two individual 
photodetectors in a small, transverse area of correlation 
equal to the coherent area of the thermal radiation than 
that of being captured in different coherent areas. For a 
large angular sized thermal source the spatial correlation 
is effectively within a physical "point" . The point-to- 
point near-field spatial correlation of thermal light has 
been utilized for reproducing nonlocal ghost images in a 
lensless configuration 0] . In a temporal HBT interferom- 
eter, the temporally, randomly radiated photons seem to 
have a two-times greater chance to be measured within a 
short time window which equals the coherence time of the 
thermal field than that of being measured in two different 
coherence time windows. For a natural thermal radiation 
source, such as the Sun, its coherence time could be as 
short as femtoseconds. It seems that the photons are 
more likely to come together and arrive at the detectors 
simultaneously if we try to detect them at some special 
space-time positions. This behavior is phenomenologi- 
cally explained as "photon bunching" [3]. What is the 
physical cause of photon bunching? We believe this is 
a reasonable question to ask. It seems that there is no 
reason to have photons created in pairs from a thermal 
source. The radiation process of thermal light is stochas- 
tic. The radiated photons should be created randomly in 
the source, rather than bunching in pairs. In fact, this 
bunching effect seems even more strange for iVth-order 
correlation measurements of thermal light, in which more 
than two photo-detections are involved. Quantum theory 
predicts that N photons have N\ times greater chance of 
being captured within the same coherent time window of 
the light field than that of being captured in N different 
coherent time window of the light field 0]. 

Recently, several papers reported high order ghost 
imaging and ghost interference experiments by using a 
pseudo-thermal source [a]. In their experiments, the au- 
thors treated the electromagnetic field from the pseudo- 
thermal source as Gaussian random variables. Since for 



Gaussian random variables any moments of order higher 
than 2 can always be expressed in terms of the first and 
second order moments, they measured high order mo- 
ments of pseudo-thermal light by measuring the first and 
second order moments with one or two CCD's, then cal- 
culated the higher order moments [fj. 

In this letter, we wish to report a three-photon tem- 
poral correlation measurement of pseudo-thermal light in 
which three photodetectors are involved in a three-fold 
joint measurement. The difference between our exper- 
iment and the experiments mentioned above is that in 
our experiment we use single photon detectors to directly 
measure the high order correlation function of pseudo- 
thermal light. We use quantum theory to describe the 
light field because it reveals the underlying physics and 
is applicable to the case when the intensity is so low that 
only a few photons are in the system. The experimen- 
tal results show that the three individual photodetectors 
D\, Z?2 and D% have a six times greater chance of being 
triggered at t\ = t% = t% than that of being triggered at 
ti =/= t2 =/= £3, apparently indicating that the photons are 
emitted in "triples" by the thermal source. 
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FIG. 1: Schematic setup of the experiment. 

The schematic experimental setup is shown in Fig. [1] 
The light source is a standard pseudo-thermal source, 
which contains a CW laser beam, a fast rotating diffus- 
ing ground glass and a focal lens (with 25.4mm focal 
length). The 632. 8nm laser beam is focused by a lens, 
onto the rotating ground glass, to a diameter of <100/im. 
The coherent laser beam is scattered by the fast rotat- 
ing ground glass to simulate a thermal field with ^0.2/is 
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coherence time. The coherence time of pseudo-thermal 
light is determined by the angular speed of the disk, the 
curvature of the focused laser beam and the transverse 
distance between the pinhole and the laser beam, the 
detailed discussion can be found in Ref.Q. Effectively, 
the rotating ground glass produces a large number of 
independent point sub-sources (3 to 5 /iin in diameter) 
with independent random phases. A pinhole with a di- 
ameter size ~ 1mm is placed 800mm from the ground 
glass to select a small portion of the radiation within its 
spatial coherence area. At this distance the size of co- 
herent area is about ~ 10mm. The pseudo-thermal light 
passes through beam-splitters BS1 and BS2. The trans- 
mitted radiation is detected by photodetector D\ , the re- 
flected radiations are detected by photodetectors D 2 and 
Z?3, respectively. To simplify the discussion, we achieved 
equal intensities in the three paths by manipulating the 
transimission-reflection coefficients of the two beamsplit- 
ters and had equal distances between the light source 
and the three photodetectors. D\, D 2 , and D 3 are fast 
avalanche photodiodes working in the photon counting 
regime. The photo-detection response time is on the or- 
der of a few hundred picoseconds, which is much shorter 
than the ~ 0.2/is coherence time of the radiation. The 
output pulses from D±, D 2 and D% are sent to a three- 
fold coincidence counting circuit which provides a three- 
photon counting histogram as a function of t\ — t 2 and 
t\ — £3, where tj, j = 1, 2, 3, is the registration time of the 
photo-detection event at D\, D 2 and D3, respectively. 

The experimentally measured and simulated 3-D third- 
order temporal correlation functions are reported in the 
upper and lower parts of Fig. [21 respectively. The sim- 
ulation is calculated based on Eq. ([5]). It is easy to see 
that (1) the measured correlation function is close to the 
simulated function, and (2) the randomly radiated "ther- 
mal photons" have much greater chance of being jointly 
detected in triples when t\ = t 2 = £3 than that of being 
detected when t\ 7^ t 2 ^ t 3 . 

From Eq. (f5]) we know that in order to claim the mea- 
sured correlation is a third order effect the contrast be- 
tween the peak and the background should be lager than 
4 to 1, corresponding to a visibility larger than 60%. 
To compare the joint counting rate at t\ = t 2 = £3 
(three-photon bunching) with the joint counting rate at 
ti ^2 7^ £3, a sliced "cross section" of the measured 
3-D histogram is illustrated in Fig. [3l The plot is a 2-D 
cross section of Fig. [2fb) sliced from the top left corner 
to the bottom right corner. The contrast between the 
maximum counting rate, which occurs at ti = t 2 = < 3 , 
and the constant background is 4.9±0.25 to 1, indicating 
a nontrivial third-order correlation function with visibil- 
ity of ~ 66%, which is greater than the 2 to 1 contrast 
(33% visibility) of HBT. This result shows that thermal 
light has a much greater chance to be bunched in triples, 
rather than in pairs. The theoretically expected contrast 
is 6 to 1 (71% visibility). In addition, the single detector 




FIG. 2: (Color Online) Measured (upper, a and b) and calcu- 
lated (lower, c and d) third-order temporal correlation func- 
tion of thermal light. The 3-D three-photon joint detec- 
tion histogram is plotted as a function of tw = ti — ts and 
£23 = t2 — ti. The simulation function is calculated from 
Eq. ([5]). In addition, the single detector counting rates for 
Di, D2, and D3 are all monitored to be constants. 
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FIG. 3: (Color Online) Crosse section of the three-photon 
coincidence counting histogram, which is sliced from the top 
left corner to the bottom right corner of Fig.fSJj. The contrast 
between the maximum counting rate and the constant back- 
ground is 4.9 ± 0.25 to 1, indicating a nontrivial third-order 
correlation function with visibility of ~ 66%. 

counting rates for D\ : D 2 , and D3 are all monitored to 
be constant. The reason for observing a lower visibility 
is due to the finite size of the detector. We expect to 
achieve higher visibility by simulating ideal point detec- 
tors. 

Phcnomcnologically, we may name the nontrivial 
three-photon correlation "three-photon bunching" be- 
cause in the measurement we find that the photons do 
have 3!-times greater chance of being jointly captured by 
D\, D 2l and D3 simultaneously than that of being jointly 
captured at different time(separated by coherent time of 
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the light field). However, does this observation really 
mean that three photons have 3!-times greater chance to 
be emitted as a bunch from the thermal source? As we 
have discussed earlier, we have difficulties following the 
photon bunching argument when N takes a large value. 

In the view of quantum interference, the observed 
three-photon correlation of thermal radiation is the re- 
sult of three-photon interference, which involves the non- 
local superposition of three-photon amplitudes, a non- 
classical entity corresponding to different yet indistin- 
guishable alternative ways of triggering a three-photon 
joint-detection event. The probability of observing a 
three-photon joint-detection event is calculated from the 
Glauber theory 



(1) 

>(i 3 )£ (+) (i 3 )^ (+) (*2)£ (+) (ii)) 



where (...) denotes an expectation value operation based 
on the quantum state of the measured electromagnetic 
field. Thermal radiation is typically in mixed states. The 
generalized density__operator for a chaotic-thermal field 
can be written as 



(2) 



where p^ is the probability that the thermal field is in 
the state 



iw>=ni 



Here we have simplified the problem to 1-D (temporal) 
with one polarization. The summation of Eq. @ includes 
all possible frequency modes u, occupation numbers n u 
for the mode u> and all possible combinations of occupa- 
tion numbers for different modes in a set of {n}. The 
quantized field operator takes the following form 



counting rate of D\, D 2 , and D 3 is calculated as 

G (3) (ii,*2,i 3 ) 
oc J duduj'duj"\f(uj)\ 2 \f(uj')\ 2 \f{u")\ 2 
1 
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[g(u, h)g(uj', t 2 )g{w", t 3 ) + g(w, ti)g{w" , t 2 )g{J t 3 



+ g(w', t x )g(w, t 2 )g(u", t 3 ) + h)g(u", h)g(w, t 3 ) 
+ g(u", h)g(u;, t 2 )g{uj', t a ) + g(w", h)g(u', t 2 )g(u>, t 3 )] 



dudu'du"\f(u)\ \f(u')\ \f(u")\ 



(4) 



Eq. (jlj is the key equation to understand the three- 
photon interference nature of the nontrivial third-order 
correlation. We notice that the probability amplitude 
has the similar form of the symmetrized wavefunction of 
three identical particles. The six terms of superposition 
in Eq. ([4]) correspond to six different yet indistinguish- 
able alternative ways for three independent photons to 
trigger a three-fold joint-detection event (see Fig. |4]). At 
t\ = t 2 = t 3 (under the condition of Z\ = z 2 — z 3 ) the 
six amplitudes are superposed constructively, and conse- 
quently G^(t\, t 2 , t 3 ) achieves its maximum value when 
summed over these constructive interferences. On the 
other hand, at t\ ^ t 2 7^ i 3 , the six amplitudes are 
not superposed constructively , the cross terms equal to 
zero and consequently G^ 3 ' (ii, t 2 , t 3 ) achieves lower val- 
ues. It is the three-photon interference that caused the 
three randomly distributed photons to have 6-times more 
chance of being captured at t± = t 2 = t 3 than that of be- 
ing captured at t\ ^ t 2 ^ t 3 . 

To simplify the calculation, taking the spectral func- 
tion f(u>) a constant within the narrow bandwidth Alu 
of the pseudo-thermal field, the normalized third-order 
correlation function g^(ti, t 2 , t 3 ) is approximately 

g {3 Hh,t 2 ,t 3 ) 



. 2r Aw(ti - 1 2 ) 
: 1 + Slnc [ _ ] 
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du f(u) g(u,tj, Zj) (v), 



(3) 



where f(to) is the spectral distribution function of the 
field, g(u>, tj,Zj) is the Green's function which propagates 
each mode of the field from the source to the jth detector. 

Assuming equal distances between the source and the 
photodetectors, i.e., z\ = z 2 = z 3 , the joint detection 



It is easy to see that when t\ = t 2 = t 3 , g^ (ti, t 2 ,t 3 ) = 6, 
the third-order correlation function achieves a maximum 
contrast of 6 to 1 (visibility ~71% ). Apparently, the 
randomly radiated independent thermal photons seem to 
have six times greater chance of being bunched in triples 
when we measure them simultaneously than that if we 
measure when ti ^= t 2 =/= t 3 . 

In fact, the quantum theory of light predicts that the 
randomly radiated independent "thermal photons" will 
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FIG. 4: Three independent photons a, b, c have six alternative 
ways of triggering a joint-detection event between detectors 
Di, D2, and D3. At equal distances from the source, the 
probability of observing a three-photon joint-detection event 
at (ti,t2,ts) is determined by the superposition of the six 
three-photon amplitudes. At ti = ti = t3 six amplitudes su- 
perpose constructively. G' 3 ' (ti, ti, ts) achieves its maximum 
value by summing over these constructive interferences. 



have AH times greater chance of achieving "AT-photon 
bunching" in an AT- fold joint-detection of N individual 
photodetectors. For large N, the contrast of the ATth- 
order correlation may achieve ~100%. The calculation is 
similar to that of the three-photon 0] 

G (K) (ii,,%) 

oc / du ... I dw"-'|/H| 2 ... |/(u/'~')| 2 



±=[ 9(u,h)...g(uj""-',t N )]\ 2 . (6) 



Eq. ([6]) indicates that the A^th-order correlation of ther- 
mal radiation is the result of A^-photon interference. It 
is easy to see that if all of the Af-photon amplitudes are 
superposed constructively at a certain experimental con- 
dition, the thermal radiation will have AI times greater 
chance to be measured under that experimental condi- 
tion. 

Although we have observed the nontrivial three-photon 
correlation with a contrast much higher than 2 to 1, and 
expect to have A^-photon correlation with a contrast of 
AH to 1, the quantum theory does not prevent having a 
constant counting rate for each of the photodetectors D\ , 
D2, and -D3, respectively. In fact, the counting rate of D\ , 
£?2j and D3 was monitored during the measurement and 
found to be constant, indicating a temporal randomly 
emitted and distributed stochastic emission process in 
the thermal source. According to the Glauber photo- 
detection theory, the counting rate of a photodetector 



is proportional to the first-order self-correlation function 
G^'(r, t), which measures the probability of observing a 
photo-detection event at a space-time coordinate (r, t). 
For the experimental setup of Fig. [TJ 



G {1 \r,t) = {E { -\r,t)E {+) {r,t)) = constant. 

In conclusion, we have observed the nontrivial 3- 
photon correlation of pseudo-thermal light. The photons 
are always randomly radiated from the thermal source. 
This experiment shows a nontrivial third-order tempo- 
ral correlation of chaotic-thermal light in which the ran- 
domly radiated thermal light is observed to have a 6- 
times greater chance of being captured by three indi- 
vidual photodetectors simultaneously than that of being 
captured by the same detectors at different times (sepa- 
rated by the coherence time of light field). The nontrivial 
A^th-order correlation of thermal light can be interpreted 
as the result of multi-photon interference, involving the 
superposition of multi-photon amplitudes, a nonclassical 
entity corresponding to different yet indistinguishable al- 
ternative ways of triggering a joint-detection event be- 
tween multiple detectors. 
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